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We study the propagating gravitational waves as a tool to probe the extra di- 
mensions. In the set-up with one compact extra dimension and non- gravitational 
physics resigning on the 4-dimensional subspace (brane) of 5-dimensional spacetime 
D , we find the Green's function describing the propagation of 5-dimensional signal 

along the brane. The Green's function has a form of the sum of contributions from 
large number of images due to the compactness of the fifth dimension. Additionally, 
a peculiar feature of the causal wave propagation in five dimensions (making a five- 
dimensional spacetime very much different from the familiar four-dimensional case) 
is that the entire region inside the past light-cone contributes to the signal at the 
observation point. The 4-dimensional propagation law is nevertheless reproduced 
at large (compared to the size of extra dimension) intervals from the source as a 
superposition of signals from large number of images. The fifth dimension however 
shows up in the form of corrections to the purely 4-dimensional picture. We find 
three interesting effects: a tail effect for a signal of finite duration, screening at 
the forefront of this signal and a frequency-dependent amplification for a periodic 
signal. We discuss implications of these effects in the gravitational wave astronomy 
and estimate the sensitivity of gravitational antenna needed for detecting the extra 
dimension. 
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1. Introduction 



Over decades since the works of Kaluza and Klein the idea of extra dimensions enjoys 
constantly growing interest of the physical community. Designed originally to serve the 
unification of the General Relativity and Maxwell electrodynamics it is now more fash- 
ionable to invoke the extra dimensions for solving the problems of hierarchy and vacuum 
energy. A recent review on old and new paths to the extra dimensions is [1]. It is cus- 
tomary nowadays to consider the Kaluza-Klein idea in the context of the brane paradigm. 
According to this picture, the non-gravitational physics, conveniently described by the 
Standard Model, is confined to a (3 + l)-dimensional subspace called brane of the (4 + n)- 
dimensional space-time with n compactified extra dimensions. Only gravitational field 
is allowed to propagate through the whole spacetime and probe the extra dimensions. 
If L is a typical size of the extra dimensions the gravitational physics on the brane is 
expected to be (4 + n) -dimensional on the scales of order L and smaller while the familiar 
4-dimensional gravitational interaction is supposed to emerge on the much larger than 
L scales. A standard way to study this is to look at the static Newton potential along 
the brane. At large distances r 3> L it takes the 4-dimensional form G^/Anr, with G4 
being the induced 4-dimensional Newton constant, while the higher dimensional physics 
manifests in exponentially small corrections, e _r//i . 

In this paper we give up staticity and look at propagation of waves. The higher- 
dimensional corrections to the purely 4-dimensional law of propagation then are no more 
exponentially small but have a power law. This makes the waves a more convenient tool 
for probing the extra dimensions. Some puzzling surprises are however awaiting for us on 
this venue. Most surprising is the fact that the causal propagation of wave signals from a 
source is radically different for odd and even dimensions. This difference* is most transpar- 
ent in the retarded potentials. In even dimensions, (3 + l)-dimensional space-time being a 
particular example, the retarded potential is proportional to the delta-function of (t — r) 
and thus is concentrated on the future light-cone and vanishing everywhere else. In con- 
trast, the retarded potential in odd spacetime dimensions has support everywhere inside 
the light-cone! Additionally, there appear non-integrable divergences on the forefront of 
the propagating wave. The above properties make the wave propagation in 5-dimensional 
space-time very much different from the propagation in the 4-dimensional case. In par- 
ticular, the night sky would look quite differently if we lived in 5-dimensional Universe. 
Even already dead stars would still produce flashing spots on the sky so that night would 
not be as dark as in our home Universe. A natural question arises^: how after all these 
differences the large scale propagation on the brane still manages to be 4-dimensional even 
though it is built out of the 5-dimensionally propagating waves. In this paper we address 
this question and moreover study the manifestations of the 5-dimensional character of the 
wave propagation in the form of corrections to the purely 4-dimensional law. 

*The mathematical theory of the propagation of waves in various dimensions is known for quite long 
time, see for example [2]. 

§ This question should have been asked years ago. However, with few exemptions [3], [4], [5] this 
property of wave propagation remains unknown in the current literature on extra dimensions. 
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We find a few interesting effects. First, we consider propagation of a signal of finite 
duration T emitted by a source on the brane. In purely 4-dimensional picture an observer 
at some large distance from the source would see a pulse of same finite duration T. In the 
5-dimensional picture with compact extra dimension the signal detected by an observer is 
a superposition of signals coming from large number of images additionally to the original 
source. In fact, the purely 4-dimensional picture arises as a result of this superposition. 
The corrections are however important. Most significantly, when the entire signal has 
passed through the point of observation there still remain signals coming from the images 
which are inside the past light-cone of the point of observation. These signals keep arriving 
and result in a tail, a power-like decaying with time signal. The whole effects thus is due 
to the combination of the 5-dimensional feature of the propagation inside the light-cone 
which we mentioned above and the compactness of the fifth dimension manifested in 
presence of large number of images. On the other hand, the shape of the forefront of the 
signal is modified by what we call screening effect: corrections appear with negative sign 
and smooth out the otherwise a step-like forefront profile. Another effect appears for the 
periodic signal. We find that the amplitude of the signal at the point of observation is 
amplified with a frequency-dependent amplification factor. Both effects are proportional 
to L in some power and are clear indications of the presence of the extra dimension. 
This makes it interesting to reconcile these effects with modern observational possibilities 
available at LIGO or LISA. 

This paper is organized as follows. In section 2 we formulate the set-up and find various 
representations for the Green's function describing propagation of the 5-dimensional signal 
along the lower-dimensional subspace (brane). We focus on the retarded potentials in 
section 3 and derive the 5-dimensional potential consistent with the boundary conditions 
we impose along the fifth dimension. The propagation of the finite pulse is considered 
in section 4 and the periodic signal in section 5. In section 6 we estimate the predicted 
effects and compare them with the available sensitivity of the gravitational detectors. 

2. Many faces of the Green's function 
2.1. The set-up 

We consider flat 5-dimensional space-time covered by coordinates X A = (x^,y) where y 
is the coordinate along the fifth dimension and x^, fi = 0, 1, 2, 3 are the four-dimensional 
coordinates. We restrict the coordinate y to lie within the interval between y = and 
y = L and we impose the Neumann boundary conditions on the perturbations of the 
five-dimensional gravitational field at the ends of this interval. L is thus a size of the 
extra dimension. Our four-dimensional world with experimental devices, standard model 
particles and human beings resign on the four-dimensional subspace at y = so that 
only gravitational field may propagate into the fifth dimension. The subspace y = (as 
well as the one at y — L) can be also thought of as a zero-tension brane. We sometimes 
will use this terminology. A generalization to branes with non-vanishing tension (the 
Neumann boundary condition then has to be replaced by a condition of the Robin type) 
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is possible although technically more involved. Another possible setting is to impose the 
periodic boundary conditions (f)(x,y = 0) = <j)(x,y = L), where (j)(x,y) stands for the 
relevant five-dimensional perturbation. In this case we have a standard Kaluza-Klein 
configuration. Although all our results are naturally extended to the case with periodic 
boundary conditions, we will stick to the scenario with the Neumann boundary conditions. 

2.2. Kaluza-Klein representation 

It should be noted that in what follows we systematically ignore the tensor structure of 
the gravitational perturbations assuming that the field equation of interest takes the form 
of the scalar type field equation 

(-d 2 + d 2 + d 2 yL ) ( f ) (y,x ^) = (2.1) 

subject to the boundary conditions 

d y (j)(y, x°, x) \ y=0 = d y (f)(y, x°, x) \ y=L = . (2.2) 

The normalized modes forming the basis of functions in the ^/-direction are 

, n = (2.3) 
nn 

cos{m n y), m n = — , n>0 (2.4) 

while in the x-direction the modes are the usual plane waves e ipx appropriately normalized 
on the delta-function. In terms of these modes the 5-dimensional propagator (or the 
Green's function) defined as a solution to the field equation with delta-like source 

-V* D 5 (X,X')=5(X,X>) (2.5) 

takes the following five- dimensional momentum space representation 

C5( X,X0 = ^/^|e^^M, (2 „ 

where p 2 = —p^ + p 2 . In Lorentzian signature there are various types of propagator 
depending on the contour chosen when the integration over momenta is taken. The par- 
ticular choices of the integration will be discussed in the next section. Here for simplicity 
of the presentation we switch to the Euclidean signature by means of the analytic contin- 
uation p = ip^ so that p 2 becomes a positive quantity and the problem of the integration 
contour does not arise. 

We see that m n introduced in (2.4) have the meaning of the Kaluza-Klein masses. 
Note that in case of periodic boundary conditions we would get a similar form for the 
propagator with different set of Kaluza-Klein masses m n = 2im/L. 

The main object of interest in the present paper is the reduction of the propagator 
(2.6) to the subspace y — 0. In other words, we are interested in the situation when the 
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gravitational 5-dimensional signal propagates between two points on the subspace where 
we happen to live. This propagation is described by the quantity 



D(x : x f ) = D 5 (X : X')\ y=yf=Q . (2.7) 

It is clear from equation (2.6) that D(x, x') reduces to the series of 4-dimensional massive 
Green's functions of the Kaluza-Klein tower 

i r oo ip(x-x') i oo 

^^^/^^"^'-^ ' (2 ' 8) 

where 

^/ 9, ,n m Ki(m\x — x'\) ,, r, rr 

D(m 2 \x-x')= - u , 1 — b-x 7 = J(x-x') 2 2.9 

47T 2 |X — X'\ v 

is the 4-dimensional propagator of massive field given in Euclidean space by the MacDon- 
ald function of the first order. Equation (2.8) gives us the first, very natural from the 
Kaluza-Klein perspective, representation for the brane-to-brane propagator as mediated 
by the set of massive Kaluza-Klein fields. 

2.3. Momentum space representation 

Another way to look at the expression (2.8) is to consider it as a momentum space 
representation for the entirely four-dimensional propagator 

D(x,x') = j^j- J d* P D(p) e**"*') (2.10) 

with the non-trivial propagator D(p) as a function of the 4-dimensional momentum p. 
D(p) is expressed in terms of the infinite sum over Kaluza-Klein masses as in (2.8). For 
the tower of these masses given by (2.4) (tracing back to the boundary conditions (2.2)) 
the sum over n can be calculated explicitly and the answer is 

D (P) = 1^ —. 2= • wr W P=VP ■ 2-11 

„_ QO p 2 + m 2 n p smh(Lp) v 

Notice that (2.11) is valid for the Euclidean signature. Switching to the Lorentzian sig- 
nature is provided by the analytic continuation p — > ip so that the Lorentzian form of the 
propagator D{p) at timelike momenta is the following 

B„M = ^f4, P=^ . (2.12) 
p sm(Lp) v 

Notice that it has poles at real values of p, Lp = irn, which correspond to the Kaluza- 
Klein modes. Additionally D M (p) has zeros at Lp = | (2n + 1), n is integer. The physical 
interpretation of the zeros of the propagator is less clear 1 '. At large values of momentum 

^ These zeros guarantee the positivity of residues of the propagator at its poles or the normal non-ghost 
nature of the corresponding Kaluza-Klein modes and also follow from the duality of the Dirichlet and 
Neumann boundary value problems in braneworld physics [6]. 
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(or, respectively for very short space-time intervals) the propagator (2.11) behaves as 
D{p) ~ the typical behavior of truly 5-dimensional propagation as it is seen from a 
lower-dimensional subspace. On the other hand, at small values of momentum (i.e. at 
large intervals on the brane) the four-dimensional behavior D(p) ~ \ is recovered. This 
guarantees us that the four- dimensional physics is reproduced at large scales even though 
the underlying physics is 5-dimensional". The cross-over between two regimes happens 
at the space-time intervals s 2 (x, x') of order L 2 . Provided we were not aware that our 
four-dimensional world is embedded into a higher-dimensional space-time we would have 
to deal with the propagator (2.12) as a non-trivial function of the 4-momentum p not 
being aware about its Kaluza-Klein origin. 

2.4. Sum over images representation 

The integration over momentum in (2.10) includes the integration over angles and the 
integration over absolute value p of the 4-momentum. Inserting the exact expression 
(2.11) into (2.10) the integration over angles can be performed explicitly and gives 

i i 7 i 

D(x) = —— jdpp I D(p)J 1 (p\x\) , (2.13) 
' ' o 

where the table integral 

7T 

J d9 sin 2 9 e iy cos 9 = - J x (y) (2. 14) 

o y 

has been used. Taking into account that on the real axis the following relation holds 
between the Bessel and Hankel functions 

Uz) = \ (ffW(z) + H[ 2 \z)) = \ (H?\z) + H?\-z)) (2.15) 

one can transform the integral to the integration over the whole real axis of p with the 
principal value prescription at p = 

oo oo 

Jdpp-D(p)J 1 (p\x\) = -P J dpp-D(p)H{ 1 \p\x\) . (2.16) 

-oo 

By closing the contour of integration in the upper-half plane of the complex variable p 
one finds that the result reduces to the contribution of poles on the imaginary semiaxis, 
p n = im n , n > (the pole at n = contributing only one half of the residue in view of 
the principal value integration above) 

OO 

f , cosh(Lp) , , ,, 1 27r ^ 

l«>«sm Jim) -m + v S-^w • (2 - 17) 

"For a curved de Sitter brane embedded in various 5-dimensional backgrounds this transition was 
analyzed in detail in [7] 
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where K\(z) = —nH < i\iz) is the modified Bessel (MacDonald) function. In view of the 
relation Ki(z) = —d z K (z) between the MacDonald functions of the zeroth and first 
order, one has 



oc 

/ 



cosh(Lp) . | . . 1 2n n ~ rr / 
dp? ■ tt/t JM^ =7TT-T2 9 ^^oH , l/r , 2.18 

Smh(Lp) L|x| ^ z=7r|x|/L 





where the series of McDonald functions can be calculated with the aid of the formula [8] 

Y K (nz) = - fc + ln— ) + — + tt V ( - 1 - —) ■ (2.19) 

t^x 2 V 2tt) 2z ^ \Vz 2 + 4n 2 n 2 2nn J 1 ; 

Putting things together this gives us yet another, coordinate space, representation for the 
propagator 

1 00 1 

D(X,X') = — V 77 rrz - 9 . (2.20) 

4tt 2 n fr / 00 ((x - x'f + 4n 2 L 2 ) 3 / 2 v ' 

This result obviously reproduces the recovery of the Green's function along the brane 
by another method - method of images. Indeed, it can be obtained by summing the 
contributions of images of the original source, located at the infinite sequence of points 

oo 

D(x) = 2 £ D 5 (X-X n )\ (2.21) 

X n = (0, 2nL), (X -X n ) 2 = x 2 + 4n 2 L 2 , (2.22) 
mediated by the five- dimensional Green's function** 

w^'w ux-xri* ■ (2 ' 23) 

The representation (2.20) of the propagator expresses the obvious fact that there many 
ways to propagate signal between two space-time points on the brane: it may propagate 
entirely along the brane (the n = term in the sum (2.20)) or through the fifth dimension 
reflecting off the second brane at y = and coming back to our brane. Each such 
propagation contributes to the sum in (2.20). Since the number of reflections can be 
arbitrary high the sum has infinite number of terms. 

When (x — x') 2 is much larger than L 2 the sum in (2.22) can be approximated by the 
integral, J2 n = I dy where y = 2Ln. Evaluating the integral we find 

1 f+°° 1 11 

D(x, x') ~ ^ J ^ dy ^ x _ x , )2 + y2)3/2 = ( x _ x ,y> ( 2 ' 24 ) 

so that 4-dimensional Green's function is recovered in this limit. It is actually a realization 
of the well-known in mathematics method called descent method to obtain Green's func- 
tion from a Green's function one dimension higher. Thus, the 4-dimensional propagator 
is recovered for large (compared to L) intervals on the brane as a result of superposition 
of contributions from infinite number of images. 

**The coefficient of two here arises from the fact that the magnitude of the source effectively doubles 
when the source approaches the brane and merges with one of its primary images - this happens with 
all pairs of higher-order images merging at X n = (0, 2Ln), n^0. 
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2.5. Newton potential 



We finish this section with a brief illustration of how the Newton potential emerges from 
our analysis. Indeed, the Newtonian potential on the brane can be easily recovered from 
(2.20) by integrating over time (it is again a realization of the descent method discussed 
in the previous subsection) 

V(r) = J dx'Dfrx') = JL £ 1 , r 2 = x 2 . (2.25) 

-oo n=-oo 

In view of the summation formula (2.11) it equals 

and, thus, reproduces at large distances r ^> L the usual four- dimensional Newtonian 
potential with Yukawa-type corrections in the four-dimensional world. For small distances 
r <C L it obviously features the five- dimensional behavior ~ 1/r 2 . 



3. Retarded potentials 

The retarded Green's function which is an object of our prime interest here can be obtained 
from the Euclidean Green's functions of the previous section by the analytic continuation. 
This can be done directly in the final result (2.20) in the coordinate representation, or 
separately for massive Green's functions of the Kaluza-Klein tower in (2.8). Here we 
choose the first method, and its equivalence to the second one is shown in the Appendix 
A. 

To distinguish the coordinates of the Lorentzian spacetime x = (x°,x) from those of 
the Euclidean space we supply the latter by the subscript E, Xe = (x 4 ,x). They are 
related by the Wick rotation 

xe = (x 4 ,x) = (ix°, x), (3.1) 
x\ = (x 4 ) 2 + x 2 = -(x ) 2 + x 2 + ie = x 2 + ie , (3.2) 

which in particular determines the analytic continuation from the (massive) Euclidean 
Green's function 

(m 2 - V 2 e)D e (x e ) = S(x E ) (3.3) 

to the Feynman propagator in Lorentzian spacetime - the boundary value on the real axis 
of the Lorentzian spacetime interval x 2 when approaching from the upper half complex 
plane of x 2 according to (3.2), e — > 0, [9] 

(m 2 - is - V 2 )D F (x) = 5(x) , D F {x) = iD E {x E ) N (3.4) 

Xe = \t% ) x ) 

(the factor i here follows from the formal identification 6(x) = i5{xe) matching with 
(3.1)). 
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As is well known [9], other Green's functions in flat spacetime can be obtained from the 
Feynman propagator by taking its real and imaginary parts and multiplying it with distri- 
butions reflecting the needed causality properties. In particular, the retarded propagator 
reads 

D R (x) = 2 6(x°)ReD F (x) . (3.5) 

We illustrate this prescription for massless field propagator. Starting with the massless 
Euclidean Green's function in four-dimensions 

D E {x E ) 2 n = (3.6) 

one finds in view of (3.3) and (3.5) the well-known retarded Green's function of a massless 
field 

^ , x 0(z°)t, i 6(t)8(r-t) n , n „, 

m 2 =o 27T 2 x 2 + ie Air r 



where in the last line the well-known formula 

1 1 



x 2 + ie x 2 — ie 



-2niS(x 2 



has been used. Important feature of this propagator in four dimensions is that it has a 
support on the past light cone of the observation point. As we shall see in a moment the 
behavior of the retarded potential in five dimensions is rather different. 

Application of this rule to 5-dimensional Green's function (2.23) or (2.20) is trickier 
because of the branching points in this expression. Wick rotation for the Euclidean 
Green's function (2.20) is based on the relation 

(x 2 ±ie) x = \x 2 \ x e ±l7TXe ^ x2) , |A|<1 . (3.8) 

Because of the restriction on A we cannot apply the above formula immediately to (2.23) 
in which A = 3/2. Therefore, we first effectively lower the degree of the singularity by 
noting that 



1 Id 1 



(3.9) 



(x 2 + ie) 3 / 2 r Or (r 2 — t 2 + ie) 1 / 2 

and then apply (3.8) to get 

( 5 ) . 1 % 1 d 9(t — r) 

(r ' t] = 4^ (r 2 - t 2 + tefl 2 = ~4^V d~r it 2 - r 2 ) 1 / 2 ' ( } 

We should apply now this formula to each term in (2.20) so that finally the kernel of the 
retarded Green's function takes the form 



. , 1*19 6(t - Vr 2 + 4n 2 L 2 ) 

D R (t,x) = -— ]T -— \ V ===>- . 3.11 
2vr i n f^ oo r Or ^t 2 - r 2 - 4n 2 L 2 
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An alternative derivation of this result, starting directly with retarded propagators for 
massive Kaluza-Klein modes in (2.8), is given in the Appendix A. Two drastic differ- 
ences from the 4-dimensional case (3.7) should be mentioned. First, the retarded Green's 
function (3.10), (3.11) has support on entire region inside the future light-cone (in 4- 
dimensional case the support was just on the surface of the light-cone) of the delta-like 
source at r = 0, t — 0. Second, the presence of non-integrable singularity on the light- 
cone in (3.10) should be noted. Additionally, all images of the original source, belonging 
to the interior of the past light cone of the observation point, contribute to the Green's 
function (3.11). 

4. Finite pulse 

The brane-to-brane Green's function (3.11) gives on the brane y = the solution of the 
equation 

V^(X) = -G 5 J(X) (4.1) 

with zero initial conditions in the past and with the source J(X) located on the same 
brane 

J(x,y)= f(t)6V>(x)6(y) . (4.2) 
Here G 5 denotes the five-dimensional gravitational coupling constant. This solution reads 

oo 

$(t,x) = G 5 / dt'D R {t-t'^)f{t') . (4.3) 

— oo 

Substituting (3.11) we have 



t—r„ 



^-^3Ml dt 'j^h • (4 - 4) 

-oo 1 > 'n rn=(r 2 + 4 n 2 i2) i/2 

When the source fit) in (4.2) has the form of a constant pulse of a finite duration T, 

f(t)=6(t)-6(t-T) , (4.5) 
the expression (4.4) (after the integration over t' has been performed) takes the form 

Ht,x) = ^^t6(t-r)I(a,P\[a])-(t^t-T) , (4.6) 
in terms of the function I(a, j3 \ [a]) of the two parameters 



a = , (47) 
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Here the square brackets in the notation [a] denote the integer part of a > and the 
function of three arguments I (a, (5\m) is given by the following finite sum 

I(a,d\m)— V , , m<a . (4.9) 

Let us first discuss the contribution of the forefront of the pulse - the first term of 
(4.5). With t reaching the value of the radial coordinate r at which the observer is 
located, the observer starts getting the signal encoded in the first term of (4.6). As long 
as this term depends on the integer part of the parameter a, m — [a], this signal is not 
continuous in time. For any value of time when a becomes integer, a new pair of terms 
at n — ±[a] appear in the sum (4.9) and the value of the sum changes discontinuously. 
This corresponds to the fact that the signals from the two new images located at y — 
±2[a]L reach the observer. Moreover, immediately at a — [a] these signals are singular, 
because they are proportional to l/yfct — [a]. This is a direct consequence of the structure 
of singularity of the retarded Green's function in five dimensions. In contrast to four 
dimensions its support is not restricted to the past light cone, rather it is given by the 
entire interior of the latter (see Eq.(3.10)), and the Green's function kernel in (3.10) has 
a singularity ~ l/(t 2 — r 2 ) 3//2 at t — > r + 0. Because of the square root the corresponding 
singularity in the m — ±[a] terms of (4.9) is not strong, and when averaged over the time 
lapse between two consequitive singularities gives, as we show below, a finite contribution. 

Immediately after t reaches the value of r only the zeroth term of the sum (4.9), 
l/a(3 2 = 8L 3 /r 2 \/t 2 — r 2 contributes to the signal (4.6). This corresponds to the contri- 
bution of only the source located at y = (and its image merging with the source and, 
thus, doubling its value - see discussion in Sect. 2). Therefore this contribution is essen- 
tially 5-dimensional by its nature and reflects the fact that the retarded potential (3.10) 
in five dimensions has non-integrable singularity on the lightcone. As time goes on more 
signals from images start to arrive so that at sufficiently large time the signal detected on 
the brane has the shape of a sequence of singular spikes. Obviously, an observer on the 
brane with devices at hand is not capable to resolve the spikes so that the time averaging 
which we employ below in the paper becomes a meaningful procedure. 

The recovery of the 4-dimensional nature of the signal for a given r takes place after 
some time when the signals from many images reach the observer and form the cumulative 
effect of many terms in the sum (4.9). This situation corresponds to the limit of a ^> 1. 
The calculation of the sum (4.9) in this limit is presented in the Appendix B and gives 
the following result in the first two orders of the asymptotic 1 / a-expansion 



vr cothvr/5 2 3 / 2 [J 

I(a,(3\ [a]) 



. S=a-[a]+l 



(3 y/a 2 + p 2 a 2 + (3 2 y a 
2 1 / 2 1 / 1 

+ 0(^1 ■ (no) 



« 2 + P 2 y/a(a - [a]) 



a 



Here 1(6) is the following integral of the parameter S — a — [a] + 1 belonging the finite 
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range 1 < 5 < 2 

This integral is a decreasing function bounded in this range by its very small value at 
5=1, 1(1) ~ 0.03 and in what follows gives a negligible contribution. 

Similarly to the situation with small a, the last term of (4.10) becomes singular also at 
large integer values of a, which corresponds to the arrival of the signal from the new distant 
pair of images located at \y\ = 2[a]L ~ t ^> L. Interestingly, though, t hat the si ngular 
spikes occurring for late times with the period At ~ 4(aL 2 /t)Aa = 2L^Jl — r 2 /t 2 — > 2L 
have now integrable singularities, so that the averaging over time of this spiky signal 
gives a finite average value. This averaging obviously implies that the temporal resolution 
threshold of the observing device is much bigger than the period 2L of spikes for hypothetic 
values of L belonging to a millimetre range. 

Thus, the temporal averaging of the spike singularity in (4.10) is finite and equals 

(we are averaging here only the most rapidly varying singular factor while the other factors 
for a 3> 1 can be regarded nearly constant). Similarly in the second term of (4.10) the 
averaging gives 

N +1 2 

'a-[a] + l)= / dayja-[a] + l = - (2 3/2 - l) . (4.13) 

[a] 

Assembling these results together and taking into account the smallness of 1(5) < 0.03 in 
(4.10), one has the dominant contribution to the temporal average of I(a,/3 \ [a]) 

7T COthTT/j 2 3 / 2 5 - 2 5 / 2 

Substituting this expression to (4.6) and taking into account the values of parameters 
(4. 7)- (4. 8) we finally have the retarded potential averaged over the time lapse between 
singular spikes 

G> ''l ' —^-) OH r) 

e irr/L _ I J V ) 
1/4 



I ( a >P\W) = ^-t^ + .™„2 + ™ ^r- » a>>1 ^ 




(2 5 / 2 - 5) 9(t - r) 

t 2 -r 2 ^>L 2 . (4.15) 



Here G 4 is an effective 4-dimensional gravitational constant 



G 4 = ^ (4.16) 
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and the subtracted term at t — T gives the contribution of the back front of the pulse (4.5) 
of finite duration T. It should be noted that the first term in (4.15) is exactly the signal 
expected in 4-dimensions. Remarkably, the radial function there is exactly the Newton 
potential which we calculated in section 2. The second term in (4.15) is a correction which 
originated entirely due to the contribution of signals arriving at a given point from large 
number of images. 

When t > T + r both signals from the forefront and the backfront (as well as from their 
images) of the pulse start to arrive. As a result, the first "static", i.e. dependent only 
on the radial coordinate r from the source, term cancels out. If we were in 4-dimensions 
this would be the only possible signal at t > T + r: after the entire pulse has passed 
through our point our devices detect no signal at all. In five dimensions the life is more 
interesting: the whole region inside the lightcone contributes to the signal detected at a 
given point. So that even after the shining candle dies off one can still see its glimmering 
light. In the present case an observer on the brane will detect a signal coming not only 
from the source but also from a large number of its images. As a result, we observe an 
interesting phenomenon: even when the pulse has passed through the point of observation 
our devices still detect some signal. This "tail" signal originates from the second term in 
(4.15) and takes the form 

(«( t ,»)) M = (** - 5)|± ( (f - T)({t L ^ Tf _ r2)1/4 - ■ (4.17) 

For later time when t>r and i>Twe have a simple expression for the tail 



(*(*,: 



=(2 5 / 2 _ 5 ) £i ^— i- . (4.18) 

tail V > 27T 2 t 5 / 2 ^ ' 



It is proportional to the square root of the size of extra dimension so that it is a rather small 
quantity. It is also a power-like decaying function of time and for late times, satisfying 
the relation t ^> r, is distance-independent. 

Finally, there is one more interesting effect: deformation of the forefront of the signal. 
It is due to the second term in equation (4.15). Noting the negative sign of this term it 
effectively leads to a relative screening of the purely 4-dimensional part of the signal. The 
typical shape of the resultant signal is shown in Fig.l. In order to analyse this effect it is 
useful to rewrite eq.(4.15) once again 

($(*,x))~$ (t,x) + $i(*,x), (4.19) 

*o(*,x) = |±, (4.20) 
, , G 4 2 5 / 2 - 5 2 2 5 / 2 - 5 ^ . , 

^ = -^- r l^m=-T,Me^ Mt ^' (421) 

where we assumed that r <t < r+T so that no contribution of the backfront is yet present, 
in terms of parameter^ a(t) = (t 2 — r 2 ) 1 ^ 2 /2L having the meaning of number of spikes 

t^We neglect the difference between a and [a] for large a. 
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Figure 1: Typical shape of the signal (for illustration we take the initial profile to have 
unit amplitude and duration T = 5, "size" of the extra dimensions is L — 0.1): forefront 
screening and the tail after backfront. 

that have already passed through the observer by the time of observation. This number 
is supposed to be sufficiently large to guarantee the formation of the four- dimensional 
signal (the first term in eq.(4.15)). It also should be relatively small for the relation t ~ r 
to hold. Eventually, the minimal value of a(t) is related to the time resolution of the 
measuring devise: it may not resolve a single spike, but should resolve the lapse between 
the arrival of the forefront and the observation moment t. Similar regime exists at the 
backfront tail: for t close to r + T the tail signal (4.17) is dominated by the number of 
spikes arriving right after the original signal ends. Later in the paper we will discuss 
possible observational implications in the gravitational astronomy of the tail effect (4.18) 
and, especially, the correction effect (4.21). 

5. Periodic signal 

In this section we consider the situation when the source on the brane radiates a periodic 
in time signal. So that we have 



where uj is the frequency of the signal, to be substituted in (4.4). The t-integral then is 
taken explicitly to produce 



/(*) = e 1 



(5.1) 




(5.2) 
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so that the gravitational field takes the form 



+00 



$(r,i) = --t-e^-dr Ho\uVr 2 + 4n 2 L 2 ) . (5.3) 

° m r n=-oo 

A possible way to estimate the infinite sum is to replace it by integral. This approximation 
is valid when u L « 1. Then we have that 

y H^\ojVr 2 + An 2 L 2 ) ~ 2 / dxH^ (uVr 2 + 4L 2 x 2 ) = — e~ iwr , (5.4) 

nt^oo Jo LU) 

where in the passing to the last expression a table integral available in [8] has been used. 
Thus, to the leading order the gravitational signal is 



G 



$(r,t)~$ (r,t) = ^ : e i ^ . (5.5) 

It is exactly the shape of the purely 4-dimensional signal arising from a periodic source. 
Thus, by replacing the sum over images in (5.3) with the integral we have explicitly 
demonstrated our earlier statement that the 4-dimensional physics is recovered from 5- 
dimensional as a superposition of signals from infinite number of images. In order to get 
a correction to the four- dimensional result (5.5) we have to estimate a deviation of the 
integral in (5.4) from the sum. First, we estimate this deviation for a single term in the 
sum (by expanding the integrand to first term of its Taylor series), 

H®(ujVr 2 + 4n 2 L 2 ) - dx H { G 2) {00^ r 2 + An 2 x 2 ) 

Jn 

= T +1 dx (H^icuVr 2 + 4n 2 L 2 ) - H$\wy/r* + AL 2 x 2 )) 

J n 

2nL 2 u ( 2 ) 



t , |r , 2 H[ 2 \u;Vr 2 + AL 2 n 2 ) . (5.6) 

Summing now this deviation for all n and again replacing the sum by the integral (which 
turns out to be exactly integrable in view of the relation H\ '(z) = {d/dz)H Q (z)) we find 

[°° dx , T L y T 2 g| 2) (Wr 2 + 4L¥)4gf 2> W • (5.7) 
Jo y/r 2 + AL 2 x 2 2 

Thus we find that the correction $i(r, t) to the purely 4-dimensional result (5.4) is 

$r(r,0 = ^:-e^^ 2) M , (5.8) 
8m r 

We are interested in a limit when the value of ur is large. Then we can use the asymptotic 
formula for the Hankel function and get 



*i(r,t) ^\ -^§^ {t - r)+il = J—e^M^t) • (5.9) 

This part of the signal is entirely due to the extra dimension. Notice the phase shift 
with respect to the 4-dimensional signal (5.5). Also, the amplitude of (5.9) now depends 
on the frequency of the signal. These two circumstances help to distinguish the purely 
4-dimensional part of the signal from that of due to presence of the extra dimension. We 
see that one of the manifestations of the fifth dimension is some frequency- dependent 
amplification of the amplitude \§ (r,t) + $i(r, t)\ of the signal from a periodic source. 
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6. Can gravitational antenna detect the extra 
dimension? 



Our analysis shows that there are two, at least in principle observable, effects of the 
compact fifth dimension on the propagation of the gravitational signal. First, it is the 
tail effect discussed in section 4: after a signal of finite duration T has passed through 
the observer there still remains a tail signal (4.18) entirely due to the presence of extra 
dimension. The late-time tail signal is a power decaying function of time, t _5//2 . In 
the four-dimensional physics the tail does not appear at all. Therefore, its observation 
would be a clear evidence for the fifth dimension. The amplitude of the purely four- 
dimensional signal at distance r from the source is the usual Newton potential (4.20), 
®o{t, r ) — G^/Aitr. Measuring the amplitude of the tail (4.18) with respect to $ we find 
that the upper bound on this ratio is given by the quantity 

L 1/2 T , , 

where we took the observation moment t > r. 

Another effect comes for a periodic signal of frequency u from a source placed at 
distance r from the point of observation. The signal is amplified due to multiple winding 
around the compact extra dimension, moreover the correction component (5.9) has a 
phase shift useful for observational purposes. The amplification factor is 1 + i] where rj 
according to (5.9) behaves as 

(6 ' 2) 



in terms of the wavelength A = 2-k/uo of the signal (in our units the speed of light c = 1). 

The two effects suggest that the extra dimension in principle could be observed in 
gravitational wave experiments and thus may have some consequences for the gravita- 
tional astronomy. The effects governed by parameters x or V are however expected to be 
extremely small. Here we give some estimates adjusted for the parameters (see for exam- 
ple [12] and [13]) of the LIGO and LISA detectors. The sensitivity of the LISA detector 
covers a low frequency range 10~ 4 — 1 Hz. For a source placed at distance of 10 Mpc from 
the Earth and producing a signal of finite duration T ~ 10 4 s (associated with the inverse 
of the lower bound of this frequency band) we have, assuming that the size of the extra 
dimension L ~ lO^ 1 cm, that x ~ 10~ 25 . The LIGO detector is sensitive for a higher 
frequency, 1 — 10 4 Hz, and is more appropriate for the observation of the amplification 
effect for a periodic signal. For the frequency lu ~ 100 Hz and the distance of 10 Mpc 
we find that rj ~ 10~ 18 . The observation of the amplification effects thus requires less 
amplitude sensitivity of the gravitational antenna. We should not forget, however, that i] 
(or x) should be further multiplied by the amplitude of the purely four- dimensional signal 
to give the value for the amplitude of the interesting part of the signal. The resultant 
effect then is too small to be experimentally observable now or in coming future. It is 
nevertheless worth noting that the late time tail (4.18) does not depend on the distance 
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from the source. Therefore, there might be a cumulative effect of superposition of the 
tails from large number of various sources. The resultant signal is proportional to the 
number of sources (both recent and in very distant past) in our part of the Universe and 
thus may have a much larger amplitude and can be singled out by the characteristic t -5 / 2 
behavior. 

There is yet another range of observation intervals which gives more hope for exper- 
imental detection. This range corresponds to (4.19)-(4.21) - the transition regime from 
the moment of the arriving forefront of the signal at t — r to late times t > r. In this 
regime the time is sufficiently larger than r to generate a large number a(t) of signal 
spikes - signals from images that already reached the observer - but not big enough to 
relax to negligible small tail rapidly decaying with time. This number a(t) should be high 
enough to guarantee the formation of the four-dimensional signal (4.19) (remember that 
it gets completely recovered at infinite t as a cumulative effect of an infinite number of 
images), but it should not be too big so that antenna would be able to resolve the lapse 
after arriving the forefront and the observation time. The correction to the 4-dimensional 
signal is given by the product a(t) -1 / 2 $o(^! x )- The lower bound on a is related to the 
time resolution At of the antenna as [min a] ~ (x) 1 ^ 2 (x) 1 ^ 2 - This gives an estimate 
oT x l 2 ~ 1CT 8 (adjusted to the frequency 10 4 Hz available at LIGO) and it makes certainly 
a much stronger effect than the ones that have been just discussed. 

Thus, the idea of echoing extra dimension(s) by the gravitational-wave astronomy 
looks feasible, and we are planning to return to this issue in subsequent publications. 

Acknowledgements 

The authors are grateful to Slava Mukhanov for his stimulating criticism. S.S. would also 
like to acknowledge a lively discussion with Dejan Stojkovic and Andrei and Valeri Frolov. 
A.O.B. is grateful for hospitality of the Physics Department of LMU, Munich, where a 
major part of this work has been done under the support of the grant SFB375. The work 
of A.O.B. was also supported by the Russian Foundation for Basic Research under the 
grant No 02-02-17054 and the grant LSS-1578.2003.2. Research of S.S. is supported by 
the grant DFG-SPP 1096. 



17 



Appendix A. Alternative derivation of the retarded 

Green's function 



The retarded Green's function of the massive Klein-Gordon operator in four dimensions 
has a representation in terms of the Bessel function of zeroth order [10, 11] 

G R (m 2 \t,x) = --l-^\j (mVW^)9(t-r)} . (A.l) 
Aixr or 1 J 

Therefore, the retarded version of the equation 2.8 involves summation over the tower 
of Kaluza-Klein masses m n = nir/L. This summation can be done with the aid of the 
Eq.(8.522) of [8] 

oo c\t m i 

E M™s/L) = — E 7=2 A2T2 > ( A - 2 ) 
n=-oo n n= _ m Vs 2 - 4n 2 L 2 

where m is the integer satisfying the inequalities 

2m < s/L < 2(m + 1) . (A.3) 
This can be rewritten for s = y/t 2 — r 2 as 



V J (nnVW^/L) = — V d{t ~ ^ + 4n2L2) 



(A.4) 



so that the resultant expression for the retarded Green's function is given by equation 
(3.11) in the main text. 



Appendix B. Calculation of /(a, (3 \ m 

The calculation of the finite sum (4.9) is based on the following representation of this sum 
in terms of the contour integral in the complex plane of the summation parameter z 

T . 7T COtll7T/3 1 f , COt7T,2 

/( °- P 1 '"' = 73 V^TW + 2iJ dz (y + 2 < W a' - ? ■ (R1) 
The closed contour of integration C m here consists of the four lines, 

c m = c+ u c m u c+°° u c-°°, 

: z = ±m + iy, — oo < y < oo, 
C ±oc : z = x±ioo, -m<x<m , (B.2) 

forming a rectangle gone around in anti-clockwise way. The contour integral in (B.l) gives 
the sum of contributions of residues at the poles of cot ttz which reproduce the needed 
finite sum (4.9) and the sum of two residues of the integrand at z — ±i/3 compensated by 
the first term on the right hand side of (B.l). 
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The integrals over C ±oc obviously give zero contribution, while the integrals over 
C+ U C~ in view of the relation 



cot 7r(±m + iy) = —i coth ny = — i ^1 + 



(B.3) 



can be reduced to the sum of the following contributions of the integral over the segment of 
the real axis, — m < z < m, the residue of the integrand at z — +i/3 and the thermal-type 
integral 



Ti 1 1 



COt TTZ 



Tt 



(p 2 + z 2 )Va 2 - z 2 



dx 



/V« 2 + P 2 J-™ (p 2 + x 2 )Va 2 -x 2 

oo 

-^J d y^r^[{f(y)-f(-y) 



(B.4) 



where the function f(y) equals 

f(y) = 

Since 

dx 



(p 2 + z 2 )Va 2 - z 2 



(B.5) 



z=m+iy 



-m ((3 2 + x 2 )Va 2 -x 2 /V« 2 + f3 2 
we finally have 

ti coth7r/9 — 1 



arctan 



V« 2 + /^ 2 



m 



(3 V< 



I (a, (3 | m) — 



2 (V^TW 

= arctan - 



cr — 



m 



(B.6) 



P V« 2 + P 2 PVa 2 + P' 2 

oo 

-^fdy^^{f(y)-f(-y)) 



P \J a 2 — m 2 



(B.7) 



So far it was exact calculation without any approximation made. We are now interested 
in the limit of large a. The calculation of the last integral in (B.7) for a ^> 1 follows from 
the asymptotic expression for the combination of functions /(±y) in (B.7) 



f(y) - f(-y) = 



^5a 2 + P 2 \ l+y 2 /5 



(l + y 2 /5 2 y/ 2 -1 



1 + 



a 



1/2 



(B.8) 



where the parameter 5 = a — m is of order unity (bearing in mind that m will be 
identified with [a]). Thus the integral in (B.7) after rescaling the integration variable 
y — > yS becomes 



oo 

-2ifdy^^{f(y)-f(-y 



5 



a 2 + P 2 \ a 



1(5) 1 + 



a 



1/2 



(B.9) 



where 1(5) is defined by Eq.(4.11) 
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In the expression for the retarded potential (4.6) the parameter m should be taken 
coinciding with [a]. Then for integer values of a the parameter 5 will be vanishing and 
1(5) will be divergent. To isolate these divergences we disentangle from the sum (4.9), 
I (a, (3 | [a]), the potentially divergent terms at n — ±[a], 

2 

I(a,f3\[a])=I(a,P\[a]-l) + . . (B.IO) 

(P 2 + [a] 2 )^a 2 -[a] 2 

The last term here becomes singular at integer values of a, but the first term is absolutely 
finite, because its corresponding value of 5 = a — [a] + 1 belongs to the range 1 < 5 < 
2. Therefore, using in (B.IO) the equation (B.7) with m — [a] — 1 and expanding the 
arctangent 



(3yja 2 + (3 2 



TT 



arctan 




one finally obtains the expression (4.10) used in Sect. 4. 



(B.11) 



References 

[1] V. A. Rubakov, "Large and infinite extra dimensions: An introduction," Phys. Usp. 
44, 871 (2001) [Usp. Fiz. Nauk 171, 913 (2001)] [arXiv:hep-ph/0104152] 

[2] G. E. Shilov, Generalized functions and partial differential equations, (Gordon and 
Breach, New York, 1968). 

[3] D. V. Galtsov, "Radiation reaction in various dimensions," Phys. Rev. D 66, 025016 
(2002) [arXiv:hep-th/0112110]. 

[4] V. Frolov, M. Snajdr and D. Stojkovic, "Interaction of a brane with a moving bulk 
black hole," arXiv:gr-qc/0304083. 

[5] V. Cardoso, O. J. C. Dias and J. P. S. Lemos, "Gravitational Radiation in D- 
dimensional Spacetimes," Phys. Rev. D 67 064026 (2003), [arXiv:hep-th/0212168]. 

[6] A. O. Barvinsky and D. V. Nesterov, "Duality of boundary value problems and 
braneworld action in curved brane models," Nucl. Phys. B 654, 225 (2003) 
[arXiv:hep-th/02 10005]. 

[7] M. K. Parikh and S. N. Solodukhin, "De Sitter brane gravity: From close-up to 
panorama," Phys. Lett. B 503, 384 (2001) [arXiv:hep-th/0012231]. 

[8] I. S. Gradstein and I. M. Ryzhik, Table of Integrals, Series, and Products, (Academic 
Press, New York and London, 1965). 



20 



[9] B. S. De Witt, Dynamical Theory of Groups and Fields, (Gordon and Breach, New 
York, 1965). 

[10] F. W. Byron and R. W. Fuller, Mathematics of classical and quantum physics. Vol. 
2. (Addison- Wesley, 1970). 

[11] R. Gregory, V. A. Rubakov, S. M. Sibiryakov, "Brane worlds: The gravity of 
escaping matter," Class. Quant. Grav. 17 4437 (2000), [arXiv:hep-th/0003109]. 

[12] B. Barish, LIGO Overview, NSF Annual Report, 23 Oct. 2002 

http:/ /www. ligo.caltech.edu/LIGO_web/conferences/nsf_reviews. html 

[13] S. A. Hughes, "Listening to the universe with gravitational- wave astronomy," Annals 
Phys. 303, 142 (2003) [arXiv:astro-ph/0210481]. 



21 



